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1 Travel time and hodograph 

Let E be an Euclidean space with the hne element ds and O C -E" be a closured 
bounded domain with C^-boundary F. Let n be a positive C^-function in the 
closure VL (refraction coefficient); consider the conformal metric g = n^ds^ 
in n. Take some points x,y E Vt and consider a geodesic curve 7 between 
these points. The length r (7) of 7 is called travel time. It is a multi-valued 
function r (x, y) of the end points x and y of 7, if, at least one geodesic has 
a pair of conjugated points. This is typical situation, if n is not constant, 
moreover, conjugated points appear for any n 7^ const, if n is constant at 
infinity, see jS] . The restriction T of the travel time function to F x F is called 
the hodograph. The branches of the hodograph T (x, y) can be distinguished 
by values of the tangent vector 9 at one end, say x. This vector can be found 
from the differential of the corresponding branch of the hodograph, see Sec. 7. 

Problem: whether the refraction coefficient n in is uniquely deter- 
mined from knowledge of its hodograph T. 

In other words, if two conformal metrics in VL have equal hodographs, 
does it imply that the metrics coincide? 

This problem was studied by Muhometov ^U], Muhometov and Romanov 
|12j . Beylkin [3], Bernstein and Gerver PPj2]- The stability of reconstruction 
of the metric from hodograph was proved in these papers under assumption 
that the travel time r is uniquely defined for any pair of points x, y G 0, that 
is the hodograph is a usual function. Then the equation ri (x, y) = T2 (x, y) 
for all pairs of points on the boundary implies rii = n2 in f2. See also the 
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papers of Croke [H], Sylvester and Uhlmann ^3]. More references are given 
in the surveys [7] and [12]. 

We prove here an estimate similar to that of fU] ('^ = 2) without assump- 
tion of uniqueness of the travel time. 

2 Conformal metric 

Take the problem in a more general setting. Let f2 be a closed bounded 
domain in = with the boundary F of class C^. Let g = Qijdx'^dx^ be a 
Riemannian metric in Q of class C^; x^jx"^ are linear coordinates in E. Let 
n gC^ (Q) be a positive function, called refraction coefficient; consider the 
conformal Riemannian metric g = ng. 

Assumptions on a refraction coefficient: (i) n G (Q) , 
(ii) for any point x G any geodesies 7 that starts at a point x & Q 
reaches the boundary F with a non tangent direction. This implies that 
there is no waveguide in Q (non-trapping geometry). 

Let T (E) be the tangent bundle on E, S (Q) be the bundle of unit circles 
in T with respect to g. Let 7 = 7 (x, 6) be the geodesic curve that starts 
from a point y = y {x,6) G F and arrives at the point x with the unit tangent 
vector 9. The parameterization 7 = 7 (x, 6*) , {x, 9) E S (Q) can be applied to 
all geodesies in Q and the function 



is defined and is C^-continuous on S (Q) . It vanishes for outgoing geodesies, 
that is for {^19) < 0, where u (x) means the inward conormal to F at x. For 
incoming geodesic we have {i^\9) > 0. So the data of travel times r (x, 9) for 
X G F is the hodograph of the metric g. 

3 Differential of the travel time 

Let n be a refraction coefficient in Q that fulfil the above assumptions (i) 
and (ii). Fix an arbitrary point {x,9) G S (Q) ; let 7 (x, 6') be the geodesic 
of the metric g = n^g that arrive at x with the (yf-unit tangent vector 9 and 
?/ G F be its initial point. Fix a vector 9 and consider r (x, 9) as a function 
of X only. 




nds 
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Proposition 1 For any point (x, 6) ^ S we have 



where e (x, 6) is a unit covector such that (e (x, ^) \6) > 0. 

Proof. Write d^r (x, ^) = An (x) e (x, 0) for a unit covector e and a 
scalar A > 0. We have {dxT (x, 6*) |^) = n (x) > according to the ((H), which 
implies A = (e (x, 6) \6)~^ . ► 

We have < (e (x, 6) \6) < 1 and can write (e (x, 0) \0) = cosu for some 
angle u G (— 7r/2, 7r/2) . 

Choose an orientation in E; it induces an orientation in the tangent bun- 
dle T {E) . Take an orthonormal basis (a, b) in the cotangent bundle T* (Q) 
that is consistent with this orientation. For a tangent vector 6 we denote 
by 6 the dual covector; its coordinate expression is 6i = gijOK Write this 
covector in the form 9 = cos (pa + sin 06 for some angle G S*^, which is a 
coordinate in the bundle of all unit covectors in Q. The coordinate system 
x^,x^,0 is globally defined in this bundle. Let d be the exterior differential 
in the complex of differential forms on S* (Q) ; write d = d^. + d^, where 
d^f = /'d0, /' = df/d(l). Differentiating yields 

dd^T = d^d^T = n (x) 

Set fj = —sin 0a + cos 06. The covectors 6,fi form a positively oriented 
orthogonal frame and satisfy the equations 

9 Afi = dV, d^ = r) A d0, 
fiAdfi = 9Ad9 = dcpdV, 

where dV^ = (det {(^jj})^^^ dx^dx^ is the Riemannian volume form in {^,g) . 
Fix the angle u by the equation e = cos a; 6' + sin u; 17. 

4 Comparing two hodographs 

Now we estimate difference between two positive functions ni,n2 (refrac- 
tion coefficients) in terms of hodographs of two conformal metrics gi = 



d^e (x, 



coscu 



smcu 



cos^ 00 



d^uj A e (x, 9) 



(3) 
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2 = n^g. Denote by 71 (x, 9) , 72 (x, 6^) the corresponding geodesies, 
by Ti{x,6), T2{x,6) the corresponding travel times. By Proposition Q we 
have dxTj{x,9) = nj (x) cos~^ ujjSj {x,9) for some angles (x'i,ix'2 such that 
= cosci;j6' + sinu jfj. Set p{x,(f)) = T2 {x,9) — Ti {x,9) and calculate the 
product 

R = dp A d^dp = (dr2 - dn) A (d0dr2 - d^dri) 
Lemma 2 We have 



R 



n2 



ni 



2 cos (a;i — UJ2) nin2 



cosci;2/ \cosuJiJ COSUJ1COSUJ2 
+ (n2 tan UJ2 — ni tan Ui) {112 — ni) A9 A fj 



9Ad9 



(4) 



Proof. Substituting 
262 A de2 



yields 



R 



n 



^ cos^ci;2 



n 



, ei A dei 

COS^ UJi 



nin2- 



62 A dei + ei A de2 



COSCUi cos 0^2 



- nin2 

We have for e = , j 



sincjidu; 



COS^ Ui cos LU2 

1,2 



sin co'2dci;2 

cos Ui COS^ UJ2 



A ei A 62. 



e A de = |^cos u9 + sin cur^j A d ycos uj9 + sin uo-f) 

cos ijj9 + sincjr/j A ^coscjd^ + smudfj — sinu;dci;6' + costudtur) 

= cos^ u;6' A d^ + sin^ ujt) Adf] + 9 A T)duj 
= 9 Ad9 + 9 Afjdu = {I + uj')9 A d9, 

since 6' A dr} = d^ A 17 = 0, and 

ei A de2 = ^coscji^ + sintuir/j A d ^cosu;26' + sinti;2r/^ 

= cos {uJi — UJ2) 9 Ad9 + cos {ui — UJ2) 9 A fjduj2 
= cos {uJi - CJ2) (1 + UJ2)9 A d9, 
ei A 62 = — sin (ui — UJ2) dV. 
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Therefore 
R 



eAde 



n 



^ COS^ UJ2 



n 



^ COS^ tUi 



nin2- 



cos {uJi — UJ2) {2 + U2 + uj[) 



nin2 sm [UJ2 — uJi) 



cos Ui cos ci;2 
cos UJ2 sin tuiu;^ — cos ui sin ct;2a;2 
cos^ a;2 cos^ Ui 



Two terms containing ^-derivatives gives 



cos {uJi — UJ2) {UJ2 + Uj[) 
COS Ui COS UJ2 



sm 1(^2 — ^1] 



cos ct;2 sin uoiuj'^ — cos cui sin 002^02 
cos^ ci;2 cos^ cji 



t^2 



cos^ cos^ UO2 ' 



times nin2, which finally yields 

2 



R 



n2 



6 Add \c0suJ2 



ni 



2 cos (cui — a;2) nin2 



COS^ 0^2 COS^ Ui 



COSCJi 

nin2 



costal COS 0^2 



UJr, 



cos^ u;i cos^ 002 



This completes the proof. 



5 The basic inequality 

Theorem 3 If ni, 112 fulfil the conditions (i) and (ii) and rii = n2 on F, the 
inequality holds 



n2 



ni 



d0dV" < 



dxP A 



(5) 



COSCi;2 COSCiJiy Js^Jt 

Proof. The last term of (j3)) is equal to the exact form 
[n2 tan (^2 — '^1 tancui) (n2 — rii) 6 At) 

since (i[d At)^ = d (a A 6) = 0. It is continuous up the boundary T as 
well as its primitive function (n2tanc(;2 — nitancui) (n2 — rii) in spite of the 
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functions tan cui, tan 1^2 tend to infinity at T. This follows from vanishing of 
111 ^ 112 on r. Therefore integration over S (Q) gives 



dp A d(f,dp 



s{n) 



Sir) 
> 



n2 



cosu;2 
n2 



ni 



cos LJi 



2 COS {uji — UJ2) nin2 



Sir) 



cosa'2 
n2 



ni 



ni 



By Stokes' 



5(-p-) \cosa;2 cosa;i 



dp A d^dp 



cos UJi 

2 

d^dV 



cos LUi COS LU2 

2nin2 

cos LJi COS LJ2 



d(l)dV 



dcpdV 



'Sin) Jsin) 
since d^d = — dd^,., and © follows. ► 



d (d^.p A d^p) 



Sir) 



d^p A d^p, 



Corollary 4 // the hodographs of rii and 112 coincide for all incoming direc- 
tions (x, 6*) , X G r, then n2 = rii. 

Proof of Corollary. It is easy to see that the condition implies 
equation rii = 112 on F. Therefore we can apply Theorem |21 and conclude 
that p{x,(j)) = 0. This yields n2Cosu;i = nicosci;2 for all x E Q and unit 
vectors 9. Fix x and choose 9 in such a way that the geodesic 71 (x, 9) has 
maximal length. The geodesic 71 is orthogonal to F at the initial point yi and 
by the Gauss Lemma (see jlj) the end points z of geodesies 71 {z, 9) such that 
Ti (z, 9) = Ti [x, 9) run over a curve / that is orthogonal to 71 (x, 9) . Therefore 
ei = ^ and oji = 0. The equation cos cji = 1 implies 112 (x) < rii (x) . The 
opposite inequality is also true. ► 



6 Application to the geodesic transform 

Corollary 5 For an arbitrary metric n that fulfils (i) and (ii) and any real 
-smooth function f in Q that vanishes at F we have 



COS^ 00 



dx9 A d^g 



(6) 
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where 

g{x,e)= [ fds. 

J-yix,e) 

Proof. We may assume that f E {VL) and apply theorem|21to rii = n, 
n2 = n + e/, where e is a small parameter. The geodesic curves for riiand \\2 
are the same up to O (e^) . Collecting the terms of order e in we obtain 
®. ► 

Example. Consider the Euclidean metric g = ds^ in the unit disc Vt. By 
a direct calculation 

J cos'^ uo \1 ~ fI / 

Write the line integrals g in the the standard parameterization G {p, (f) = 
g [x, 9) of line integrals, where 

X = (cos s, sin s) , —7i/2 < s < 7i/2,6 = (cos (f), sin </)),s — 7r/2<</)<s + 7r/2, 

p = cos {lp — (f)) , (f = (p + IT /2 

and 

g', = - v^t^^g;, g'^ = v^t^g; + g; 

The inequality © takes the form 

y (l^^j f (^)dx<- ds g'Ax,e)g'^{x,e)d<j> 

= j^^ (v^T^ (g;)' + g;g;) dpdy.. 

7 Calculation of angular derivative of the travel 
time 

Take an arbitrary refraction coefficient n in f2 that satisfies (i) and (ii) and 
determine the angular parameter ^ of a ray 7 in terms of the hodograph 
function. Let s be the arc length along F that increases clockwise. 
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Proposition 6 Let (x, 6') , x G F be an arbitrary outward direction and 7 (s) 
be the family of geodesic that joins some point y G F and x (s) G F with 
outgoing tangent vector 6 (s) defined for s G [0,e] (or for s G [—£,0]^ such 
that X (0) = X, 6 (0) = 6. Then we have 

simp — — — \s=Q, 

n [X) us 

where ip is the angle between the outward conormal —v (x) and the outward 
covector 9. 

Proof. By the Gauss Lemma the tangent vector 9 is orthogonal to the 
circle of points with the g-distance r (x, y) from y. ► 

This formula also holds for any point x conjugate to ?/ in 7 (0) . In par- 
ticular, 9 = —u, if dr/ds = 0, that is the geodesic 7(0) arrives to x with 
the normal direction —u (x). The angle of ^ of the ray 7 (x, 9) is equal to 
(p = ip + P + 71, where (3 = arg u. 
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